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The two-loop electroweak radiative corrections to the muon's anomalous magnetic moment, = (g^ — 2)/2, 
are presented. We obtain an overall 22.6% reduction in the electroweak contribution from 195 x 10~^^ to 
151(4) X lO"^'^. Implications for the full standard model prediction and an upcoming high precision measurement 
of are briefly discussed. Some aspects of the calculations are discussed in detail. 



O 



> : 

cn . 
o\ ■ 
m ■ 

^ ; 

o 

\o 

5 ■ 
i!^ ■ 
Oh: 

I 

Oh; 
^ ■ 



X 



The anomalous magnetic moment of the muon, 
= ((7^ — 2)/2, provides both a sensitive quan- 
tum loop test of the standard SU (3)c x SU (2) l x 
U{1) model and a window to potential "new 
physics" effects. The current experimental av- 
erage 11 



116592300(840) x 10" 



(1) 



is in good agreement with theoretical expecta- 
tions and already constrains physics beyond the 
standard model such as supersymmetry and su- 
pergravity |2,^, dynamical or loop muon mass 
generation compositeness leptoquarks |^ 
etc. 

An upcoming experiment E821 
at Brookhaven National Laboratory is expected 
to start taking data in 1997. With one month of 
dedicated running, it is expected to reduce the 
uncertainty in a'j^P to roughly ±40 x 10^^^, more 
than a factor of 20 improvement. With subse- 
quent longer dedicated runs it could statistically 
approach the anticipated systematic uncertainty 
ofabout ±10-20x10"" g. At those levels, both 
electroweak one and two loop effects become im- 
portant and "new physics" at the multi-TeV scale 
is probed. Indeed, generic muon mass generating 
mechanisms (via perturbative or dynamical loops 
[Q) lead to Aa^ « m^/A^, where A is the scale of 
"new physics". At ±40 x 10"^^ sensitivity, A « 5 
TeV is being explored. 

To fully exploit the anticipated experimental 
improvement, the standard model prediction for 
must be known with comparable precision. 
That requires detailed studies of very high or- 
der QED loops, hadronic effects, and electroweak 
contributions through two loop order. The con- 
tributions to Qa are traditionally divided into 
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^QED ^ 116584706(2) x 10"". 



(3) 



The uncertainty of the QED result is well within 
the ±20- 40 X 10"" goal. 

Hadronic vacuum polarization corrections to 
enter at 0(a/7r)^. They can be evaluated via a 
dispersion relation using e^e" — > hadrons data 
and perturbative QCD (for the very high energy 
regime). Employing a recent analysis of e^e~ 
data ||l^,^ along with an estimate of the leading 
0{a/^ effects, we find Q 
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= (vac.pol.) = 6934(153) x 10" 
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Unfortunately, the error has not yet reached the 
desired level of precision. Ongoing improve- 
ments in e'^e" — > hadrons measurements at low 
energies along with additional theoretical input 
should significantly lower the uncertainty in (^. 
Nevertheless, reducing the hadronic error below 
±20 X 10"^^ or even ±40 x 10"^^ remains a 
formidable challenge. 

The result in (^) must be supplemented by 
hadronic light by light amplitudes (which are of 
three loop origin) pO|-p^. Here, we employ a 
recently updated study by Hayakawa, Kinoshita, 
and Sanda O] which gives 
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(light by light) = -52(18) x 10" 
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However, we note that the result is somewhat 
dependent on the low energy model of hadronic 
physics employed and continues to be scrutinized. 
Combining (^ and (||) leads to the total hadronic 
contribution 
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QED loops have been computed to very high or- 



Now we come to the electroweak contributions 
to , the main focus of our work and the impetus 



for forthcoming experimental effort. At the one 
loop level, the standard model predicts Mi 



a^W(lloop) 
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where = 1.16639(1) x 10"^ GeV-^, M = Mw 
or Mniggs, and the weak mixing angle sin^ 9w = 
s"^ = l-M^/Mj = 0.224. We can safely neglect 
the O {ml/M'^) terms in (0). 

The one loop result in (j^) is about five 
to ten times the anticipated experimental er- 
ror. Naively, one might expect higher order (2 
loop) electroweak contributions to be of relative 
0{a/7r) and hence negligible; however, that is not 
the case. Kukhto, Kuraev, Schiller, and Silagadze 
(KKSS) have shown that some two loop elec- 
troweak contributions can be quite large and must 
be included in any serious theoretical estimate 
of or future confrontation with experiment. 
Given the KKSS observation, a detailed evalua- 
tion of the two loop electroweak contributions to 
is clearly warranted. We have reported the 
complete results of such an analysis in two re- 
cent papers |^,^. The types of diagrams con- 
tributing to two-loop electroweak corrections are 
shown in ^ These fall into a number of differ- 
ent categories. Diagrams with a closed fermion 
loop, together with a class of bosonic diagrams, 
represent corrections to vertices and propagators 
of the one-loop electroweak diagrams (shown in 
the first line in Fig.|l|). 

In addition we have to calculate new types of 
diagrams which appear at the two-loop level, such 
as mixing, an induced jjH vertex etc. (lines 
two and three of Fig.|]). Finally, there are non- 
planar diagrams and diagrams with quartic cou- 
plings; some examples are shown in Fig.^. 

The final result of the total electroweak contri- 
butions to the g-2 of leptons - including one- and 
two-loop corrections - is 

„EW _ on in-- 14 



,EW 



= 3.0-10' 
= 151-10"" 
= 4.7-10""^. 
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In this paper we discuss some details of the pro- 
cedures we have employed in that project. In par- 
ticular we discuss the method for projecting out 
the Pauli formfactor in d dimensions, and some 
aspects of the fermionic contributions. 




Figure 1. Bosonic diagrams contributing to g-2 
of the muon 





Figure 2. Examples of non-planar diagrams and 
diagrams with quartic couplings. 



Diagrams with charged bosons 

Contributions of diagrams with a fermion loop 
connected to the muon line via two charged 
bosons are shown in Fig. 3. Fermions with isospin 
±1/2 arc denoted by u and d, respectively. Of 
course, in the case of leptons diagrams (a) and 
(f) vanish. 

We first consider the case when the fermions 
in the loop belong to the first two generations. 
Here the masses of the fermions in the loop do 
not influence the result very much and we neglect 
them. The ratio of the neglected terms to the 
result is at most of the order of 
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< 0.3% 
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with nic denoting the mass of the charm quark 
1.5GeV. In the massless approximation we 
need to consider only the first four diagrams in 
Fig. 3. In ref. |26| it has been argued that the di- 
agrams 1(a) and 1(b) vanish by virtue of Furry 's 
theorem. We find that this is not true even after 
adding contributions of all fermions in a genera- 
tion. Furry's theorem consists in the observation 
that the sum of contributions of diagrams with 
two different orientations of the fermion loop van- 
ishes. This is not the case for diagrams 1(a) and 
1(b) because for every fermion flavor interacting 
with the external photon there is only one pos- 



sible orientation of the fermion line. Only those 
parts of the expressions which contain a single 75 
cancel out after adding contributions from the up- 
type quark, down-type quark and from the lep- 
ton. 

Let us consider diagrams 1(a) and 1(b) in more 
detail. The contributions of these diagrams can 
be divided up into those with an odd number of 
vector couplings to the fermion loop {VVV or 
T/j4A) and with two vector, one axial vector cou- 
pling (AVV). The dominant contribution is ob- 
tained by taking massless fermions in the loop. 
For diagrams of type 1(a) and 1(b) we find (the 
subscript indicates the fermion to which the pho- 
ton couples) 
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with 



C = 



8V27r3 
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Here we used the shorthand notation VVV for 
the sum of VVV+VAA. Ar?^'^'-'^^ = 3 is the num- 
ber of colors (with 7V^^*°" = 1 understood) and 
Qf denotes the fermion charge. Only the par- 
ity conserving part contains divergences; they are 
proportional to the lowest order ^WW vertex. 

Summing the AVV part over one light genera- 
tion (l.g.) we obtain 

^a^r = (12) 
due to the fact that in the standard model 

E^cQ/ = 0. (13) 
/ 

The VVV contribution of one light generation is 
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equal (even equal to zero as in our example) just 
like in the QED case. The Furry theorem is not 
applicable because, unlike in QED, the charge 
flow through the fermion loop destroys the sym- 
metry between the fermion loop diagram and the 
diagram with reversed fermion direction. 

We summarize now the total contributions 
from a single light generation (where appropri- 
ate, we always include a mirror counterpart of 
the diagrams) 
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These expressions have to be multiplied by a 
color factor where necessary. We put the relevant 
Kobayashi-Maskawa matrix elements equal to 1. 
Adding contributions of all fermions we obtain for 
a single light generation 



C 10 
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We now proceed to the contribution of the third 
generation. For the r lepton loop we can still use 
eq. (fBI). We obtain 



In diagrams with quarks we neglect the mass of 
the b. To leading order in — ^ we obtain for the 
sum of all diagrams containing top and bottom 
quark loops 

(18) 
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The singular terms — are cancelled by renor- 
malization of the W boson mass present in the 
one loop electroweak contributions to muon g-2. 



This example demonstrates explicitly that the 
contributions from three vector couplings to the 
fermion loop do not cancel. We would like to 
emphasize that this has nothing to do with the 
masses inside the fermion loop. We can put them 
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Figure 3: Diagrams with a charged boson and 
a fermion loop contributing to the muon g-2. 
Crossed circles denote interactions with an exter- 
nal photon. 

Projector for g-2 form factor 

An efficient projector for g-2 calculations has been 
obtained in ref.|^^. Here we extend their formu- 
lae to d = 4 — 2e dimensions. 




Figure 4: Definition of momentum flow 

We consider the most general matrix element of 
a current between spin 1/2 fermions 



{af\J^,{x)\a 

Uf{p2 



(19) 



^iW7m - ^F2{t)a^,A'' + -F3(i)A^ 
2m m 



■Ui{pi)e 



with A = pi — p2 and A^ = t. For on-shell exter- 
nal fermions we have 

p1^pI = rr? (20) 
Introducing p — \ {p\ +P2) we obtain in addition: 



p^ = -(4m^-t), 



p- A = 0. 



(21) 



Conservation of the electromagnetic current re- 
quires F-^if) = 0. Fi{t) is the charge form factor, 



F2{t) the magnetic moment form factor, G2{t) 
the electric moment form factor. The anomalous 
magnetic moment of the fermion a is given by 



(22) 



In order to extract the magnetic moment form 
factor, one can introduce a projection operator 
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In order to determine the coefficients gi we take 



the trace of TV^M'': 



Tr(iV^M^) = { [8to2 + 2t{d - 2)] gi 
+2{t - Am^)g2} F,{t) 
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This leads us to the following system of equa- 
tions for F2{t): 

= [8m^ + 2t{d-2)]gi + 2{t-Am^)g2 

it — 4to^) 

1 = 2t{d-l)gi+t^-^-^g2 (25) 

and we can put 53 = already at this step. The 
solution is 

— 2m^ 
= (d-2)t(t-4m2) 
2m2 [4to2 + {d- 2)t] 
{d-2)t{t-Am?Y 
and we find 

-2to2 
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m) = 
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It is also possible to extract directly the 
anomaly i^2(0). As a first step the general am- 
plitude can then be expanded to first order 



in A^: 



Mm(p,A) 



M^(p,0)+ A,^M^(p,A) 

^v,{p) + a,t;{p). 



A=0 

(29) 



The next step is to average over the spatial di- 
rections of A with the formulas 



(A^A,) = 



(30) 



after that, the hmit i — > can be taken. The 
resuh is 



1 



(31) 
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